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Summary 

LF and MF ground-wave propagation over a realistic terrain has 
been found to be very different to that over an assumed smooth earth. This 
report describes the basic nature of these differences, which include a non- 
simple relationship for field variation with distance and a non-free-space 
relationship between electric and magnetic fields. The theory of why these 
effects occur is summarized, and the measurements taken to demonstrate 
their nature are given. Also a description is given of a computer calculation 
method which has been devised to make use of the new knowledge. 
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Fig. 1 - A mapping of the groundwave attenuation on an Argand diagram of the parameter 
p = —jwr} 2 r/k. The values of the contours are dB relative to a perfectly conducting plane 
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List of principal symbols 



a 
e 
A 

J 

E 

S 

h 
u 

m 

Zo 
n 

R 
V b 

F 
P 



H 



Xj,x 2 ,x 3 



surface impedance relative to the in- 
trinsic impedance of free space 
conductivity (siemens/metre) 
dielectric constant 
wavelength (metres) 
relative surface impedance for a 
smooth uncluttered earth 
current density (amps/sq. metre) 
electric field strength (volts/metre) 
effective complex conductivity (sie- 
mens/metre) 

height of clutter (metres) 
dummy variable used in integration 
(see Equation (4)) 

number of unipoles per square metre 
(see Equation (5)) 
characteristic impedance (ohms) 
number of trees per square metre (see 
Equation (6)) 

resistance per metre (ohms/metre) 
relative impedance of a perfectly con- 
ducting surface containing hemis- 
pherical bosses (see Equation (7)) 
fraction of area filled by a hemis- 
pherical boss 

a parameter used to calculate surface 
wave attenuation (see Fig. 1) p = 
—}irq 2 r/X 

distance from transmitter to receiver 
(metres) 

magnetic field strength (amps/metre) 
percentage of ground covered by 
buildings 

percentage of ground covered by trees 
constants in Equation (12) deter- 
mined by empirical means 



1, Introduction 

This report summarizes results of investigations 
into ground-wave propagation in a realistic terrain. 
The work was carried out in the l.f. and m.f. broad- 
cast bands, i.e. 150 kHz to 285 kHz and 525 kHz 
to 1605 kHz. The investigation began because of 
poor coverage by a few transmitters in urban areas, 
but the scope widened when it was found that 
propagation in a realistic terrain differs signifi- 
cantly in important respects from the conven- 
tional theoretical model. 



A report 1 and an IEE paper 2 described the 
effects of a built-up area on the magnetic field 
component from an m.f. broadcast station, and 
electric/magnetic field ratios were discussed in 
Electronics Letters. 3 References 1 and 2 give 
the results of a measurement campaign in London, 
showing a surprising field strength distance 
relationship. The explanation of this was shown 
to be that metallic structures considerably modi- 
fied the surface impedance to a value which could 
not be obtained on a smooth earth. The theory 
worked out for this indicated that in clutter the 
mean value of the ratio of electric to magnetic 
field would not be as in free space. This was later 
proved by measurement and published in 
Reference 3. 

The recent work has increased the amount of 
measured evidence to include a greater variety of 
areas and frequencies. In particular the effect 
of wooded areas has been studied. Some refine- 
ments have been made to the theory, and to the 
method of calculation which allows for the 
realistic terrain. 

A receiver, aerials and analysing equipment 
were made for the measurements. These enabled 
continuous recordings and statistics of fields to 
be taken on the move in a car. In addition it was 
possible to make quick changes of the received 
frequency and the aerials (from electric to mag- 
netic detectors). Other reports will be published 
to describe this equipment. 



2. Theory 

The field at a point is a function of the relative 
surface impedance values, tq, existing between 
transmitter and receiver. For a smooth earth, 
conductivity o (siemens/metre), dielectric constant 
e and wavelength A (metres) the surface impedance 
effective at grazing incidence is 



*?.=■ 



[e-j60o\-l]* 

e-j60aX 



(1) 



Values of a and e are given by Terman 4 for 
different types of terrain. With the relative 
magnitudes quoted it is sufficient for our case 
to take 
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The error of this approximation is numerically 
insignificant compared with that due to the assum- 
ption of a smooth earth. 

In a realistic terrain, containing buildings, 
trees, surface irregularities etc., a vertical current 
density, J, will be induced by the electric field, E. 
Thus we can think of an effective complex conduc- 
tivity, vertically, being given by: 



I 

' E 



(3) 



At the top of the clutter, height h, a value for 
surface impedance was derived in Reference 1. 
This may be written as: 



This work gives results of the form : 



% = 



■FW 



(7) 



where h is the radius of the bosses (or height 
above the plane). F is the fraction of area 
filled and W is a complex number which has 
a real part at least 6 times larger than the 
imaginary and is a function of the ratio of 
skin depth to radius of boss. 

The analogy to undulating terrain is weak In two 
respects: the plane is not perfectly conducting 
but has surface Impedance ?7 S and hills unlike 
bosses have a height much less than their base 
radius. However, the analysis does suggest a 
reasonable answer might be given by: 



V = V S + 



jlirh 



du 



1 -j6GAS(u) 



(4) 



where the dummy variable u ranges from the base 
of the irregularities to the top. S is written as a 
function of u because it may vary with height. 
This leaves the problem of giving values to S in 
the practical cases. For this purpose we made 
three categories. 



a) 



c) 



Highly conducting earthed unlpoles, of less 
than quarter of a wavelength high, were 
given the values 



irmh 2 

S = j u 



(5) 



where Z is the characteristic impedance 5 of 
the unipole and m is the number of them per 
square metre. 

b) For resistive objects such as trees: 



n 
5 = — 

R 



(6) 



where R is the resistance per metre and n is 
the number per square metre. 

For terrain undulations it Is useful to con- 
sider an analysis by Wait 6,7 of propagation 
over hemispherical bosses of finite conduc- 
tivity sitting on a perfectly conducting plane. 



2tt 



(8) 



where y is a real number approximately equal to Vi. 

The result obtained in Equation (8) is in 
keeping with the other equations for clutter 
and with other types of reasoning. It is equivalent 
to saying that: 



5 = 



60A 



(9) 



where C is approximately equal to one. 



In many cases the term added to the smooth 
earth surface impedance, because of the cluttered 
and Irregular nature of the terrain, is sufficient 
to radically change the field levels. In fact the 
whole form of the field strength/distance relation- 
ship may be altered as can be seen from Fig. 1. 
Here use is made of the fact that surface wave 
attenuation can be calculated via a complex 
variable we call p. This variable can be formed 
from the basic parameters of surface impedance, 
distance and frequency as follows: 



p = — J7TT7 2 r/X 



(10) 



The plot is then a mapping of the magnitude 
of attenuations on the complex plane p. Thus 
a straight line from the origin corresponds to 
a radial line from the transmitter for homogeneous 
ground constants; distances on the diagram are 
then proportional to distances from the trans- 
mitter (r). For a smooth earth the relevant part 
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of the diagram is always just below the positive 
real axis. However, with the surface impedance 
of a realistic terrain many results will be obtained 
from the upper region of the plot. Here changes 
in field value are more drastic. 

As well as obtaining a result for surface im- 
pedance, the theory developed in Reference 1 
gave a result for the electric to magnetic field 



ratios in the realistic terrain. 
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(ID 



For a propagation path that changes signifi- 
cantly along its length a method of calculation 
is required which readily allows for changes of 



BROOKMANS PARK TX 




ROUTE TAKEN FOR 
DEMONSTRATION 

CHARTS 



XTO HASTINGS 



(RA-182) 



Fig. 2 - Sketch map of the London Area 
-3 - 
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Fig. 3 - Sketch map of the Birmingham Area 
-4- 



surface impedance. This may be achieved by an 
integral equation method of field calculation, also 
described in Reference 1 . 

With this theory it is possible to attempt field 
calculations for a realistic terrain. However, 
certain quantities are still rather nebulous and it 
is left to the following empirical work to put these 
on a firmer basis. 



3. The measurements 

The measurement campaigns were organized to 
cover the frequency range and include a variety of 
terrain, with particular emphasis on that terrain 
where the most dramatic effects might be 
expected. The frequencies measured were 200, 
526.5, 692, 908, 1151, 1214 and 1457 kHz. The 
areas selected for measurement were London, 
Birmingham, the New Forest, a radial from Droit- 
wich to Hastings (on the South-East coast) and a 
small circular tour in Surrey. In each case both 
the electric and magnetic fields were measured. 
Fig. 2 is a sketch map of the London area showing 



two radials from the Brookman's Park transmitter 
along which many measurements were concen- 
trated. Also shown is part of the radial from the 
Droitwich transmitter to Hastings and a route 
taken to produce demonstration charts. Fig. 3 is 
a sketch map of the Birmingham area showing 
three radials from the BBC transmitter (Sutton 
Coldfield) and another three from the indepen- 
dent local radio transmitter called BRMB (Langley 
Mill), measurements being as close to the radials 
as roads permitted. Fig. 4 is a sketch map of the 
New Forest area showing the route taken when 
measuring the Hartley transmitter. Table 1 shows 
the way in which the 'clutter' is marked on these 
maps. 

The way in which to present the measured 
and calculated results gave us some difficulty. 
The points which went into our decision were: • 

1 . We wished to draw attention to the difference 
from the free space ratio which existed in 
the ratio of electric to magnetic field. (Con- 
ventional measuring methods added a con- 
fusion here because receivers usually express 
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Fig. 4 - Sketch map of the New Forest Area 
- 5 - 



TABLE 1 
Topographical categories 



Type of Area 


Class 


%of 

ground 

cover 


Height (m) 
of clutter 


Shading on 
maps 


City centres 


A. 


30 


20 




WllllIIIIi 


Town centres 


B 


20 


15 




WMM§M 




Suburban 


C 


15 


10 




i ! 


Semi-Rural 


D 


10 


10 








Wooded 


F 


10 


15 




m^m 




Bare 


Blank 















results in juV/m even though they detect 
magnetic field.) 

2. We wished to emphasize the attenuation of 
field relative to that which would be obtained 
over a perfectly conducting plane earth. (To 
include the perfectly conducting plane-earth 
attenuation in our plots would make the 
particular points of concern less clear.) 

3. We could not accurately find the absolute 
value of the field over a perfectly conducting 
plane earth, from a given transmitter, unless 
parameters such as aerial efficiency could be 
determined. 

These points led us to choose a presentation 
which means that on the plots described below, 
electric field is expressed as 20 log 10 (Er) mV and 
magnetic field as 20 log (HOirHr)^ where r 
is the distance from transmitter to receiver. 

Fig. 5 is a plot of measurements from the 
Brookman's Park transmitter. The crosses and 
squares are median values for 500 metre stretches, 
for magnetic and electric fields respectively. The 
results are not like those obtained over a smooth 
earth, even if changes of conductivity were in- 
cluded. They do, however, follow a trend that 
would be obtained if, from Fig. 1, a line were con- 
sidered which went from the origin and passed 
near the 'hole'. 

The electric and magnetic fields are both 
plotted in terms of millivolts per metre, assuming 
the free space ratio to obtain a direct comparison. 
It can be seen that the electric field has a strong 



tendency to be a lower value, this tendency in- 
creasing with the amount of clutter. 

Fig. 6 is a plot of measurements obtained 
in Birmingham from the BBC transmitter (1457 
kHz) and the BRMB transmitter (1151 kHz). 
The measurements taken from the BBC site show 
similar trends to London results, but from the 
BRMB site there is not such a dramatic change 
with distance. This is attributed to the lower 
frequency and to the fact that the BRMB site is 
more suited to serve the town centre. 

Fig. 7 is a plot of measurements from the 
Bartley transmitter (692 kHz) taken through the 
New Forest. These were taken to check the 
theory about the tree effects; unfortunately the 
New Forest is not sufficiently dense and extensive 
to give very dramatic results. 

Fig. 8 is a plot of measurements from the 
Droitwich transmitter (200 kHz) along a radial 
line ending in Hastings on the South-East coast. 
The main points to be noted here are that the 
magnetic field decays little from the perfectly 
conducting earth value and the electric field is 
generally lower than this magnetic field value. 
The smooth earth calculation would give greater 
attenuation of the magnetic field with distance 
and does not predict the relatively reduced electric 
field. 

The theory indicates that there is no fre- 
quency factor in the relative depression of the 
electric field due to vertical metallic structures. 
To test this theory, 200 kHz from the Droitwich 
transmitter was also measured on the routes from 
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F?§-. 5 - Measurements through London from the Brookmans Park transmitter 
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Fig. 6 - Measurements taken in the Birmingham Area 
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Fig. 7 - Measurements taken through the New Forest 
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Fig. 8 - Measurements taken from Droitwich to the South-East coast 
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Fig. 9 - Comparing the relative reduction of electric field on 200 and 1457 kHz 



the Brookman's Park transmitter through London. 
The plot of Fig. 9 was then made; in this the 
crosses result from the analysis of measured results 
and a regression line (solid) was compared to a 45 
line (dashed). It can be seen that the theory 
and measurements are in reasonable agreement. 

Fig. 10 gives examples of pen recordings 
taken on the route in Surrey (see Fig. 2). The 
figure attempts to show, on short pieces of chart, 
as many features as possible but not extreme 
cases (the town included is not heavily built up 
with tall buildings). On Fig. 10(a) the route 
first crosses a motorway intersection, then des- 
cends a steep hill with a high footbridge crossing 
near the top. The town centre is at the foot 
of the hill. The sharp drop-out of electric field 
occurs when passing through a narrow section of 
street. Fig. 10{b) was recorded in an area that 
is sparsely built-up, with an uphill route and 
other features as marked. The high spikes of 
field on the 1457 kHz electric-field trace are 
associated with overhead power lines. The 
wooded area has the greatest effect on the 200 
kHz electric field, which again supports the theory 
i.e. vertical resistive objects have a frequency 
dependence such that lower frequencies have 



the greater attenuation. 

It was suspected that the attenuation of 
electric field caused by trees would be less in 
winter because of the reduction in sap resulting 
in an increase in resistance. Therefore, the tour 
was redone in winter and the hypothesis con- 
firmed. In particular, the 200 kHz section of 
chart marked 'heavy tree cover' showed an increase 
of the electric field of about 2.5 dB. Of course, 
this could be attributed to the lack of leaves, but 
the theory indicated that the total tree resistance 
is the important factor. 



4. The calculation system 

The measurements obtained have supported 
the theory put forward, but they must do more 
than that because it is not practicable to obtain 
ail the required quantities from theory. A further 
practical point is that even though some para- 
meters may be reasonably estimated it may be very 
time consuming and costly to do so. For a calcu- 
lation system to be of value it must either do 
something that cannot be done by measurement 
or it must do it quicker, cheaper and without 
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requiring excessive amounts of skill and knowledge 
from the user. With these thoughts in mind the 
topographical categories of Table 1 were devised 
and assigned as shown in Figs. 2 to 4. This 
reduces to a minimum the amount of data required 
from the user, about clutter. 

For each of a series of V2 km squares along 
the radial line we give the '% of ground cover' by 
buildings the symbol P b and by trees P t . Then 
using the information contained in Equations (5) 
and (6) we get: 

60X5(u)^X(x 1 P t +x 2 P b )+jx 3 P b u (12) 

where x , x 2 and x are to be determined by 
empirical means. 

The value of x was determined by actual 
measurement of the resistance of trees and found 
to be approximately 12 x 10 -5 . The values 
of x 2 and x were determined so that the London 
measurements were calculated with a minimum 
of error. This resulted in x 2 being given the value 



-2 



165 x 10 8 and x equalling 2.5 x 10 



For the determination of electric field via 
Equation (11) it is not sufficiently accurate to 
use the percentage of ground cover given in Table 
1. These figures are averages over large areas, but 
the electric field is modified in a way which is 
very dependent on the immediate surroundings. 
Thus if a measurement is taken in parkland the 
ratio of electric to magnetic fields will be close 
to that of free space. However, in a highly built- 
up street the electric field is reduced to a greater 
degree than the figures of Table 1 would give. 
Our measurements were taken driving down such 
streets and it was found that by doubling the 
ground cover values in Table 1 fairly accurate 
predictions of these electric fields were obtained. 
This extra problem associated with electric field 
will make its accurate prediction more difficult 
than that of predicting the magnetic field. 

The data for a calculation was then as shown 
in Fig. 11, together with explanations of the 
various items of data required. The data for 
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Fig. 11 - Data for input to Computer Program 



height of terrain is of little importance but is 
included for the sake of completeness. 

An example of the print-out of results is 
shown in Fig. 12. The output gives the wave- 
length in metres, labelled as 'LAMBDA(M)'. The 
columns are as follows. 

1. The distance of field point from trans- 
mitter in metres. 

2. The terrain height in metres. 

3. The average height assumed for clutter in 
metres. 

4. The percentage of area covered by build- 
ings. 

5. The percentage of area covered by trees. 

6. The ground conductivity in millisiemens 
per metre. 

7 & 8. The surface impedance, rj, including the 

clutter factor as defined in Equation (4). 

9. The magnitude of the loss of the magnetic 

field relative to a perfectly conducting 

flat surface. 



10. 
11. 
12. 



13. 



The phase of the loss factor. 

Column 9 expressed as dB. 

20 log 10 (1207rH x distance) mV . The 

reasons for choosing this representation 

are given in Section 3 above. 

This gives a figure for the additional loss 

(dB) of the electric field relative to the 

magnetic, i.e. 20 log (1207rH/£). In free 

space this value would be zero. 



The calculations of magnetic field are plotted 
in Figs. 13—16 together with the measurements 
so that the degree of success may be assessed. In 
order to make a comparison between calculations 
and measurements it was necessary to assume a 
field value which would be obtained on a perfectly 
conducting earth from each transmitter. This was 
done by consideration of the transmitter charac- 
teristics and the field values close to the trans- 
mitter. 

Fig. 13 shows the results from Brookman's 
Park on radials through London. The field trends 
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Fig. 13 - Comparison of measurements and predictions in London 
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Fig. 14 - Comparison of measurements and predictions in Birmingham 
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Fig. 15 - Comparison of measurements and predictions in the New Forest 
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/%. 1<5 - Comparison of measurements and predictions from Droitwich 
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are reproduced with reasonable accuracy, much 
improved relative to results without the new 
clutter factors. However, we have sacrificed the 
possibility of greater accuracy for two reasons: 

1. The simple clutter categories make prepar- 
ation of data for the programs easier, cheaper 
and quicker. 

2. The magnitude of the clutter term is restric- 
ted to keep results in a domain clockwise of 
a line, on Fig. 1, which is 51° from the origin. 
In the extreme examples of the high fre- 
quency end of the MF band and highly built- 
up areas larger angles can exist, but these 
cases are very sensitive to data accuracy 
and pass through a point of zero field. 

The effect of our constraints is particularly notice- 
able on Fig. 13(c), but we consider the smoother 
calculated result to be preferable to a dramatic 
calculation in a case where it is not realistic. 

Fig. 14 shows the Birmingham results, which 
were not used in the development of the predic- 
tion system. Again, the field trends are well 
reproduced. However, on Figs. 14(a)— (c) the 
calculations have a reduction at about 5 km, 
associated with the trees in Sutton Park, which 



is greater than measured. This may be explained 
if these trees are less tall and dense than the values 
assumed in our category. 

Fig. 15 shows the results from the Bartley 
transmitter indicating that a reasonable predic- 
tion is made for the effect of trees. In particular 
it may be noted that more attenuation of the 
magnetic field occurs on Fig. 15(b) beyond the 
10 km range where there are more trees. Fig. 16 
shows the results for the 200 kHz measurements 
from Droitwich to Hastings. In this case only 
the terrain undulation effect was allowed for in 
the calculation because it was considered to be 
the dominating influence on the field variation. 

Table 2 gives some statistics of the differences 
between calculated and measured attenuations. 
There are some rather large errors in the electric 
fields which are due to particularly dense and 
tali buildings. We cannot expect accurate calcu- 
lations in these cases when only using the simpli- 
fied clutter data. To give an idea of the advantage 
of the system proposed we have calculated the 
same statistics for the Brookman's Park 1457 
kHz transmission, on 169°, without using the 
new clutter factors. These figures are shown 
in brackets on the table. It can be seen the 
errors are much larger; this is despite the fact 



TABLE 2 





Statistics of differences 


between calculated and measured attenuations 
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that we have allowed reduced conductivity values 
for the built-up areas. 



5. Conclusions 

The work described in Reference 1 has been 
extended to include trees and terrain undulations, 
both of which have been found to have a relevant 
modifying effect on propagation of ground waves. 
Many more measurements have been made to 
substantiate the theory and produce factors needed 
for calculations. Electric and magnetic fields 
have been measured and it has been shown that 
there may be large differences from the free 
space ratio. 

The calculation system has been designed 
so that a prediction can be made easily; the 
necessity to choose heights and percentage of 
ground coverage is removed by having simple 
categories. The accuracy of the calculation is 
such that predictions can be made which are much 
better relative to the assumption of a smooth 
earth. 
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